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1 .  Introduction 


It  is  well  known  that  dual -measurement  techniques  are  needed  for 
accurate  estimation  of  rainfall  rate  using  radar.1  Of  the  several  dual' 
measurement  schemes  postulated  in  the  literature,  the  technique  based 
on  differential  reflectivity,  or  ZDR,  appears  most  promising. 2* 3»H» 5 
This  final  report  deals  with  the  theoretical  and  experimental  statis¬ 
tical  properties  of  Z^R.  The  latter  are  based  on  the  time  series  radar 
data  obtained  with  the  high  resolution,  "fast"  polarization  switching 
Chilbolton  radar  system  located  in  Southern  England  and  operated  by 
the  Rutherford  and  Appleton  Laboratories. 6 

The  major  results  of  the  research  are: 

A.  The  optimum  estimator  of  Z^R  was  derived  and  excellent  agree¬ 
ment  was  obtained  between  theoretically  and  experimentally 
obtained  standard  errors.7  Using  pulse-to-pulse  polariza¬ 
tion  switching  with  a  prf  of  610  Hz,  Z^R  defined  as 

10  log  (Zjj/Zy)  can  be  measured  with  standard  errors  in  the 
range  0.1  -  0.2dB  at  a  single  range  gate  using  about  40  -  60 
independent  samples.  This  high  accuracy  is  achieved  due  to 
the  high  degree  of  cross-correlation  between  Z^  and  Zy, 
estimated  to  have  median  values  of  0.97  or  higher. 

B.  The  effects  of  noise  on  the  accuracy  of  ZQR  measurements 
was  examined.  A  signal-to-noise  ratio  of  20dB  or  greater 
is  required  to  achieve  high  accuracy  estimates  of  ZnD. 

C.  The  time  series  radar  data  from  England  was  used  to  deter¬ 
mine  the  fluctuation  spectrum  (periodogram)  of  Z^R  for  the 
purpose  of  searching  for  signatures  characteristic  of 
polarization  dependent  backscattering  processes  such  as 
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drop  vibration,  canting,  etc.  A  theoretical  model  of 
raindrops  vibrating  at  their  natural  oscillation  frequency 
(with  no  motion)  was  developed  and  the  periodogram  of  ZDR 
was  computed  for  comparison  with  measured  data. 

D.  Earlier,  the  authors  had  suggested  using  the  differential 
propagation  phase  shift  in  rain  between  the  horizontal  and 
vertical  polarizations  as  another  parameter  to  be  combined 
with  ZpR  for  accurate  rainfall  measurement.0  The  accuracy 
of  this  differential  phase  shift  measurement  was  estimated 
as  a  function  of  the  cross-correlation  parameter  and  the 
number  of  independent  samples  in  a  manner  similar  to  that 
used  for  the  Z^R  analysis. 

This  final  report  summarizes  the  results  (outlined  above)  obtained 
during  the  contract  period.  All  the  tasks  set  forth  in  the  original 
work  statement  have  been  satisfactorily  completed. 


2.  Statistical  Model 

A  scattering  model  will  be  considered  for  which  the  co-polar 
backscattered  signal,  for  a  horizontally  or  vertically  polarized 
transmitted  pulse,  can  be  written  in  the  form 

Aj^exp  (iIh(V)  "  22  Fj(H,V)  exp  (l0j(H,V))  ^ 

where  A^  v  and  v  are  the  resultant  amplitude  and  phase,  and  Fj  ^  ^ 

and  9 . ...  are  the  backscattered  amplitudes  and  phases  from  the  jth 
J  * ) 

scatterer  in  the  scattering  (pulse)  volume.  The  subscripts  (H,V) 
refer  to  horizontal  or  vertical  polarization,  respectively. 

We  resolve  v  into  their  phasor  components  (X^  YH)  and 
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(Xy,  Yy)  where  A^y  *  (X*  y  Y2  y)**  and  assume  that  the  set  of  random 
variables  (X^f  Y^,  Xy,  Yy)  are  normally  distributed  with  zero  mean  and 
covariance  matrix  given  by: 


VH 

0 

P°HaV 

0 


P^Hav 


0 

P°HaV 


p¥v 


[2] 


where  E[X^y]  =  E[Yj|^y]  *  ^H,V  =  E(Y^Yy)  =  po^Cy. 


2.1.  Optimum  Estimator  of  Zqr 

This  model  is  valid  for  coincident  sampling  of  the  pulse  volume  at 
the  two  polarizations  (H,  V) .  and  Zp^(dB)  are  defined  as  aH2/0y2 

and  10  log  (a^/ay2),  respectively.  The  optimum  estimator  of  ZD^  is  of 
the  form 


where  Aft.^  form  a  pair-wise  set  of  identically  distributed  (Eq.  [2]) 
random  samples.  This  estimator  is  termed  the  square  law  estimator  and 
the  calculated  standard  error  (dB)  as  a  function  of  sample  size,  m,  and 
cross-correlation  coefficient  p  is  shown  in  Fig.  1.  Note  the  advantage 
*  of  having  p  as  high  as  possible  in  order  to  reduce  the  standard  error 

of  Zpj^(dB)  to  0.1  -  0.2dB  with  40  -  60  independent  samples. 


2.2.  Effects  of  Noise 

The  effects  of  white,  Gaussian  noise  can  be  simulated  by  altering 
the  covariance  matrix  as  defined  in  Eq.  [2]  to  F  ♦  !  0n2  where  I  is  the 


STANDARD  DEVIATION  (OBI 


Figure.  1.  Standard  deviation  (dB)  of  the  square  law 
estimator  as  a  function  of  sample  size  m 
and  cross-correlation  coefficient  p. 


identity  matrix  and  an2  is  the  noise  power,  assumed  equal  at  the  two 
polarizations.  The  effect  of  noise  results  in  a  reduction  of  p  to  p* 
defined  as 
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where  the  signal-to-noise  ratio,  SNR  =»  a^/a^1  and  ZDR  =  oH2/av2.  It 
is  estimated  from  Eq.  [4]  that  an  SNR  >  20d8  is  required  to  maintain 
the  high  accuracy  in  the  measurement  of  ZpR. 


2.3.  Differential  Phase  Shift 

The  potential  use  of  the  differential  phase  shift  in  rain  due  to 
propagation  depends  to  a  large  degree  on  the  possible  measurement 
accuracy  that  can  be  achieved  for  this  parameter.  The  probability 
density  function  of  A  »  -  <fry  is  given  by9 


(1  -  P2)  t(l  -P2cos2A)1j  pcosAlir  -  cos-1  (ocosA)|  ] 
f(A)  - - rr - —  [5] 

2tt(1-pWA)*  11 


where  4^  v  is  defined  by  Eq.  [1]  and  the  covariance  matrix  Z  is  defined 
in  Eq.  [2].  We  assume  phase  difference  measurements  A  are  available 
at  two  range  gates  separated  by  a  distance  l  so  that  the  mean  propaga¬ 
tion  differential  phase  shift  is  given  by  Re(k^  -  ky)2£,  where  y  is 
the  propagation  constant  at  the  two  polarizations.  Neglecting  the 
differential  phase  shift  due  to  backscatter  at  the  long  wavelength  (say 
10  cm) ,  the  fractional  standard  deviation  (FSD)  of  the  differential 
phase  shift  was  determined  to  be 
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„„„  standard  deviation  _  1  ,/T  \  f*  .2*,^  ».P 

FSD  =  - 5^ - -  2lRe(kH--kv)  Vi  (4  A2f(A)dA| 

[6] 

In  Fig.  2,  we  show  plots  of  m  versus  p  for  a  constant  FSD  of  0.1  and 
0.2  assuming  a  rainfall  rate  of  100  mm  hr*1  over  a  path  length  of 
21  =  2km.  At  lower  rainfall  rates,  the  number  of  samples,  m,  increases 
very  rapidly  if  FSD*s  of  0.1  or  0.2  are  to  be  maintained.  The  effects 
of  noise  can  be  determined  by  replacing  p  by  p*  as  defined  in  Eq.  [4]. 

3.  Radar  Measurements 

The  measurements  used  to  obtain  the  statistics  of  ZDR  were  made 
using  a  10  cm  wavelength  radar  on  a  fully-steerable  25  m  diameter 
antenna  in  Southern  England.  The  essential  characteristics  of  the  radar 
system  are  given  in  Table  1.  The  antenna  was  fed  from  its  prime  focus 
by  a  circular  scalar  feed,  so  as  to  best  match  the  aperture  of  the  dish 
and  present  the  same  beam  pattern  for  vertical  and  horizontal  polariza¬ 
tions.  A  fast  polarization  switch  coupled  this  feed  to  the  radar  in 
such  a  way  that  alternate  (transmitted)  radar  pulses  were  vertically 
and  horizontally  polarized.  The  received  power  was  always  copolar. 
Because  of  the  high  radar  power  and  the  rapid  switching  rate,  a  mechan¬ 
ical  rotating  chopping  disc  was  used  rather  than  an  electrical  switch. 
The  radar  PRF  (600  Hz)  was  generated  from  the  rotating  disc  (6000  rpm) 
to  ensure  synchronization. 

The  received  radar  power  for  each  of  the  polarizations  was  then 
processed  through  the  same  receiver,  IF  log  amplifier,  and  analog  to 
digital  converter.  Measurements  were  made  only  with  the  antenna 


Number  of  Independent  Sample, 


Figure  2 


Number  of  independent  samples,  m,  versus  cross-correlation 
coefficient  p  to  determine  mean  differential  phase  shifts 
with  Fractional  Standard  Deviations  (FSD)  of  0.1  and  0.2. 
Rainfall  rate  is  100  mm  hr"1,  and  wavelength  of  10  cm  is 
assumed. 
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TABLE  1 

RADAR  CHARACTERISTICS 


Frequency 

3076  MHz 

Peak  power 

S00  kw 

Pulse  width 

0.5  ys 

Pulse  repetition  frequency 

610  Hz 

Polarization  sampling 

alternate  vertical 
and  horizontal 

Antenna  measured  gain 

53  dB 

Antenna  3dB  beamwidth 

16  arcminutes 

IF  bandwidth 

10  MHz 

Digitization  range 

180  steps  of  0.3dB 

Range  gate  width 

0.3  us 
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pointing  fixed,  and  with  a  single  range  gate.  Two-hundred  8-bit  samples 
from  each  polarization  were  stored  together  on  magnetic  tape  as  a  data 
record,  corresponding  to  a  two-channel  stream  lasting  640  ms  (see  Fig.  3) 
860  ms  were  required  between  records  for  data  transfer  and  reset  condi¬ 
tions.  Some  3,500  records  of  this  type  were  available  for  analysis. 

Some  of  the  data  records  were  obtained  for  radar  pulse  volume  return 
from  near  the  center  of  raincells,  others  from  the  edge  of  raincells, 
and  yet  others  from  regions  of  the  melting  layer. 


4.  Data  Analysis 

The  first  step  in  the  analysis  of  the  data  was  to  conduct  a  goodness 

ru  Ah/oh 

of-fit  test  to  the  distribution  F(u)  3  JQ  f(x)dx,  where  u  3  Xyfdy  ’ 
and  the  probability  density  function  f(u)  is  given  as7 


f(u) 


2u(l  +  u2)(l  -  p*) 
[(1  +  u2)2  -  4p2u2 ] 


[7] 


Note  that  the  above  function  is  based  on  a  single,  coincident  sample  of 
A(j  y  from  a  given  range  gate.  The  Chi  Square  test  was  performed  to  see 
if  the  measured  data  originated  from  the  null  distribution,  F0(u). 
Approximately  1,500  records  were  analyzed  and  about  90%  were  accepted 
at  a  1%  significance  level.  This  high  acceptance  rate  established  both 
the  high  quality  of  the  radar  data  and  the  adequacy  of  the  statistical 
model. 

The  data  was  next  analyzed  to  estimate  the  standard  deviation,  o, 
of  2QR(dB)  using  the  optimum  square  law  estimator  defined  in  Eq.  [3]. 
The  median  value  of  a  based  oi\125  data  sets  was  found  to  be  0.214dB 
which  compares  very  well  with  theoretical  estimates  in  the  range  0.16  - 
0 . 32dB ,  using  p  in  the  range  0.92  -  0.98.  Further  range  or  time 


Periodogram  of  20  log  (Ah) 
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Figure  3.  One  record  of  time  series  data  from  a  fixed  range 
gate  lasting  640  msec. 


4800 


3600 


1200 


1200 


25 

Frequency,  Hz 


Figure  4.  Periodogram  of  20  log  (Ah)  with  adjusted  mean  based 
on  200  samples  and  a  time  duration  of  640  msec. 


averaging  is  expected  to  reduce  o  to  around  O.OSdB.  The  sane  125  data 
sets  yielded  a  median  value  of  p  *  0.97.  Additional  details  of  the 
preceding  data  analysis  are  given  by  Bringi  et  al.7 

The  measured  data  is  based  on  a  time  lag  of  1.6msec  between  the 
horizontally  and  vertically  polarized  signals,  while  the  theory  assumes 
coincident  sampling  at  the  two  polarizations.  It  is  expected  that 
coincident  sampling  would  produce  higher  cross-correlation,  thus 
further  reducing  the  standard  error  in  the  measurement  of  ZpR.  However, 
system  complexity  would  increase  due  to  a  two-channel  receiving  system 
with  more  stringent  antenna  and  feed  requirements.  Decorrelation 
between  AH  and  Ay  for  sequential  sampling  can  arise  due  to  two  causes, 
viz.,  (a)  drop  vibration,  canting,  etc.  and  (b)  the  Doppler  spectrum. 

In  severe,  convective  storms  the  decorrelation  due  to  the  Doppler 
spectrum  would  dominate  causing  a  significant  reduction  in  p(~0.90), 
while  for  a  narrow  Doppler  spectrum  the  decorrelation  would  be  slight 
due  to  drop  vibration,  canting,  etc.  As  a  result,  more  "independent" 
pairwise  samples  of  ZH  y  will  occur  due  to  the  combination  of  the  above 
effects.  Nevertheless,  sequential  polarization  sampling  is  preferred 
due  to  its  simplicity  and  adaptability  of  existing  single  polarization 
radar  systems  to  possible  future  modification. 

5.  Spectrum  Calculations 

The  time  series  radar  data  obtained  from  measurements  was  used  to 
estimate  the  fluctuating  spectrum  of  20  log  AH  and  20  log  (A^/Ay) .  A 
typical  periodogram  of  20  log  A^  is  shown  in  Fig.  4  with  adjusted  mean 
using  a  standard  procedure  called  SPECTRA.10  Almost  all  of  the  data 
records  show  a  similar  narrow  spectrum  with  estimated  standard  deviation 


of  about  0.5  msec-1  at  10  cm  wavelength. 

Sample  periodograms  of  20  log  (AH/AV)  are  shown  in  Fig.  5  and  6 
with  adjusted  means.  The  spectrum  of  Fig.  5  is  typical  of  records 
which  have  very  low  mean  ZDR  s  0  . 25dB) .  These  spectra  are  broad  as 
compared  to  those  having  higher  mean  ZpR's  (*1  -  1.5dB)  of  which  a 
sample  is  shown  in  Fig.  6.  Unfortunately,  no  data  having  mean  ZpR 
values  in  the  range  2  -  5dB,  as  typically  found  in  rain  in  the  conti¬ 
nental  United  States,  were  available.  It,  therefore,  would  be  desirable 
to  analyze  data  from  a  variety  of  rain  conditions  and  at  a  larger  number 
of  range  gates  before  drawing  any  firm  conclusions  regarding  the  meteor¬ 
ological  significance  of  the  ZpR  spectrum. 

5.1.  Fluctuation  Model 

It  is  well  known  that  raindrops  vibrate  at  their  natural  frequency 
while  falling  at  their  terminal  velocity.  This  vibration  causes  small 
fluctuations  of  ZpR  about  a  mean  value  and  it  is  of  interest  to  calculate 
the  fluctuation  spectrum  based  on  a  simple  model  with  no  motion.  An 
exponential  drop  size  distribution  was  assumed  with  each  drop  oscillating 
at  a  natural  frequency  given  by  f  *  4.22a0"1*1*7  where  a0  is  the  radius 
of  the  equivolumic  spherical  drop.  The  drops  were  assumed  to  oscillate 
sinusoidally  between  spherical  and  oblate  spheroidal  shapes  with  maximum 
amplitudes  of  20%  of  the  equivolumic  spherical  diameter.  Gans'  scattering 
theory  was  used  to  calculate  ZpR(t).  The  resultant  periodogram  with 
adjusted  mean  is  shown  in  Fig.  7  for  a  mean  ZQR  of  1.3dB.  Even  though 
the  fluctuations  of  ZpR  about  its  mean  value  is  very  small,  the  compari¬ 
son  between  Fig.  7  and  the  actual  periodogram  of  ZQR  as  shown  in  Fig.  6 
is  striking.  If  the  Doppler  spectrum  has  a  large  variance,  the  ZpR 


Periodogram  of  20  log  (Ah/Av)  with  adjusted  mean  based  on  200  samples 
and  a  time  duration  of  640  msec. 


and  a  time  duration  of  640  msec. 


raindrop  vibration  model. 
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spectrum  is  not  expected  to  yield  any  information  on  raindrop  vibration. 
Other  methods,  such  as  broad  band  noise  techniques  or  coherent  polariza¬ 
tion  radar  techniques,  may  be  required  to  obtain  further  information 
about  the  hydrometeors. 1 l» 12  Also,  higher  resolution  recording  of  the 
co-polar  reflectivities  is  desirable  in  order  to  retain  all  of  the 
important  spectral  characteristics  of  ZpR. 

6.  Conclusions 

The  differential  reflectivity  parameter,  ZDR,  shows  great  promise 
as  a  valuable  parameter  in  radar  meteorology.  It  can  be  measured 
accurately  and  rapidly,  as  shown  by  the  statistical  characteristics  of 
the  Zp£  signal  when  compared  with  actual  radar  data  obtained  from  the 
Chilbolton  radar  system.  The  optimum  estimator  of  Z^R  was  described, 
and  the  sequential  method  of  polarization  sampling  is  recommended, 
provided  sampling  occurs  within  the  decorrelation  time  of  the  scatterers 
in  the  pulse  volume.  It  is  shown  that  ZDR  can  be  measured  with  a  stan¬ 
dard  error  of  0.2dB  at  a  single  range  gate  using  40  -  60  independent 
samples  in  a  time  duration  of  1  sec.  Furthermore,  the  radar  data  is 
shown  to  be  consistent  with  a  statistical  model  describing  the  back- 
scattering  from  the  pulse  volume.  The  effects  of  noise  in  the  model 
were  also  considered,  and  the  same  model  was  used  to  predict  the 
accuracy  to  which  differential  propagation  phase  shifts  could  be 
measured. 

The  time  series  data  obtained  from  the  Chilbolton  radar  was 
analyzed  spectrally  for  both  the  horizontal  amplitude  returns  and  the 
ZpR  signal.  Large  differences  in  the  shape  of  the  spectrum  were  found 
for  conditions  of  small  (~0.2SdB)  and  large  (l.SdB)  mean  ZDR's.  The 
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fluctuation  spectrum  of  the  horizontally  polarized  signal  was  very 
narrow  with  an  estimated  Doppler  standard  deviation  of  0.5  msec"1  at 
10  cm  wavelength  for  all  the  measured  data.  The  ZpR  spectrum  obtained 
with  real  data  was  compared  with  the  expected  ZpR  spectrum  from  a  model, 
exponential  rain  drop  size  distribution  where  the  drops  are  fixed  in 
space  and  vibrate  at  their  natural  oscillation  frequency.  Although 
the  ZQR  fluctuations  in  the  model  are  very  small,  the  shapes  of  the 
model  spectra  and  the  real  spectra  (for  mean  ZpR  *  1.5dB)  were  similar. 
It  is  recommended  that  further  time  series  data  be  analyzed  under 
conditions  of  larger  mean  ZpR  and  at  a  larger  number  of  range  gates. 
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